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Abstract 



We present simple and predictive scenarios for neutrino masses which lead 
to a bimaximal pattern of mixings that apparently provides the best fit to the 
current solar and atmospheric neutrino data. It is shown that an approxi- 
mate global Le — L^ — Lr symmetry, broken only by Planck scale effects, can 
naturally explain the bimaximal mixing pattern. In one scenario, the solar 
neutrino oscillation parameters are induced entirely through radiative renor- 
malization group effects. A one parameter model describing all neutrino data 
is presented, which predicts C/gS — \/rn^/rn^ ~ 0.017. A second two param- 
eter scheme predicts Ues to be vanishingly small. We also show how these 
scenarios can be extended to accommodate the LSND observations within a 
3-|-l neutrino oscillation scheme by including a mirror sector of particles and 
forces. 

I. INTRODUCTION 



There appear to be persuasive arguments in favor of the possibility that both the solar 
and the atmospheric neutrino anomalies are explained by a specific pattern of neutrino 
mixings where both the z/g — z/^ as well as the — z/^ mixing angles are nearly maximal, 
while the — z/^ mixing angle is small • The z/g — z/^ mixing explains the solar neutrino 
deficit whereas the z/^ — Ur mixing accounts successfully for the atmospheric neutrino data. 
This pattern has been called in the literature as bimaximal mixing |^ . The best fit values for 
the mass difference squares for the atmospheric and solar data seem to be Am^^^ ~ 3 x 10~^ 
eV2 and Am|, ~ (10"^ - 10"^) eV^. 

A major theoretical challenge is to provide an understanding of this pattern in reasonable 
extensions of the standard model (SM). There are many models in the literature which 
provide ways to arrive at the bimaximal pattern starting with different assumptions. In this 



note, we suggest several simple scenarios which are quite predictive and therefore testable 
in the near ftiture in ongoing and proposed neutrino oscillation experiments. 

Before proceeding to present our schemes, we wish to motivate some special forms of the 
light neutrino mass matrix that lead to bimaximal mixings. These matrices are obtained 
under some restricted but well motivated set of assumptions, mindful of possible approximate 
symmetries that might be present in the Icpton sector. Once identified, controlled breaking 
of these approximate symmetries can lead to very predictive schemes for neutrino mixings, 
as we shall demonstrate. 

We work in a basis where the charged lepton mass matrix is diagonal. We define the 
neutrino mixing matrix (the MNS matrix) U as: 



(1) 



where the states on the left and the right hand sides denote flavor and mass eigenstates. 
The light neutrino mass matrix Aiu is given by the equation 

= UMdU^ (2) 

where Aid is diagonal and has the form = Diag {mi, 1712,1713). For the mass eigenvalues, 
there are essentially three different possibilities of physical interest: (a) hierarchical form 
where mi -C m2 <^ m^; (b) inverted form where \mi\ ~ \m2\ S> 1 777.3 1 and (c) degenerate 
form where mi = m; m2 = m + 62', m^ = m + S3, with ^2,3 <^ tti. Since experiments 
indicate that the i^fj, — Ur mixing angle is nearly maximal, and that Ue3 is very small, we 
choose the following approximate form for the mixing matrix U : 
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Here s = sin^ is the solar neutrino oscillation angle, which is of order one, while the 
corresponding angle for atmospheric neutrino oscillation has been set to 1/ In the exact 
bimaximal limit, we will have s = c = l/\/2- 

Consider case (a) in the exact bimaximal limit. Let us denote 7773 = m, m2 — S and set 
777i = 0. The matrix A4u is then give by 



/ 6/2 -5/{2V2) -6/{2V2) \ 

-6/{2V2) m/2 + 5/4 5/4 - m/2 . (4) 
\-6/{2V2) S/4-m/2 6/A + m/2j 



For case (c), the form is the same as Eq. (4), except that there is an additional piece 
proportional to the unit matrix, Diag{l, 1, 1)7770, with m^ setting the overall scale for the 
degenerate mass. The relative hierarchy among the parameters in this case will be 7770 ^ 
777 ^ 5. 

For case (b) with the inverted spectrum, we get (choosing for simplicity the case c—s — 

and 7773 = 0, 777i = 777 = —7772 + S), 
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/ 5/2 {2m- 5)/{2y/2) {2m- 5)/{2y/2)\ 

M^!/^ = {2m - 6)/{2y/2) 6/4 S/4 . (5) 

\{2m-6)/{2V2) 6/4 6/4 J 



Of these forms, Eq. (5) is particularly interesting since if we set 5 = it reduces to 



with a redefined m. The mass matrix given in Eq. (6) has the interesting property that it 
has a global symmetry Lg — — Lr {Li is the ith lepton number) along with a permutation 
symmetry 5*2 acting on the z/^, Ur flavors. This new symmetry unraveled in the leptonic 
world is perhaps analogous to the local B — L symmetry that is revealed by the seesaw 
mechanism |^ for neutrino masses. 

As it stands, A^^J'-* of Eq. (6) cannot lead to solar neutrino oscillations since Ami, = 0. 
One needs to break the L^ — L^ — Lr symmetry in order to make the model realistic. Since 



the mass parameter m in Eq. (6) is fixed by the atmospheric data i.e., m ~ yAm^j^, the 
strength of the Le — L^ — L^- breaking should be of order Amg/Am^j^ ^ (3x 10^^ — 3 x 10^^). 
The smallness of this symmetry breaking parameter suggests that Le — L^ — Lr as well as 
5*2 symmetries are indeed approximately conserved. 

In this paper, we attempt to seek realistic gauge models that realize Eq. (6) as the 
leading approximation of the light neutrino mass matrix. The correction terms to the mass 
matrix will involve very few parameters, so that there is a certain predictive power for 
neutrino mass difference squares as well as mixings. The hope is that such models can then 
be tested in future experiments and may shed light into the nature of new physics operating 
in the neutrino sector and perhaps more generally within quark lepton unified schemes. We 
present several examples of such models. First we show that global L^. — L^ — Lr symmetry, 
broken only by Planck scale effects, can give a rather good fit to the neutrino data. Then 
we present a scheme where the neutrino sector has only one parameter which is fixed by 
the atmospheric neutrino data. Radiative renormalization effects from the charged lepton 
sector then determines the rest of the parameters in the neutrino mass matrix, making the 
model extremely predictive. For instance, in this model, we obtain the relation 



where r is a calculable parameter which depends only on the scale of new physics (e.g. 
seesaw scale) governing the smallness of m^. In this model, the mixing parameter [/es is 
predicted to be t/gs — — 0.017. We also present a variation of this model with one 

more parameter, where f/es can be larger (even close to the present experimental upper 
limit of 0.2). We present a second scenario which has two parameters. This model predicts 
bimaximal neutrino mixings with f/gs = 0. Thus any evidence for nonzero f/es would rule 
out this scenario. Finally, we consider extensions of these schemes that may provide a 
simultaneous fit to solar, atmospheric as well as the LSND neutrino oscillation data. As is 
well known, this needs the introduction of a sterile neutrino. We achieve this in a manner 
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that introduces only two more parameters into the theory by invoking the mirror matter 
models, where all new parameters (except the weak scale and possibly the QCD scale) are 
fixed by an exact mirror symmetry 



U. LE-Lf,- Lr SYMMETRY AND BIMAXIMAL MIXING 

Suppose that Le — — Lr symmetry is a global symmetry of the Lagrangian broken 
only by Planck scale induced terms. (Quantum gravity is expected to break all non gauge 
symmetries.) The Lagrangian in this case will have the general form: 

£ = i-(aLe0L^0 + 6Le0i:.0) + ^(L,0L^0) + /i.c. (8) 

IVl IVl PI 

where are fiavor indices. Here is the standard model Higgs doublet. M is the seesaw 
scale, which can be near the GUT scale ~ 10^^ GeV, so the first set of terms in Eq. (8) will 
dominate over the last set. Note that the Lg — — symmetry is broken by the Planck 
scale induced terms. The leading contribution to Aiiy is of the form in Eq. (6), with the 
two nonzero entries unequal in general (denote them mi and m2). Upon diagonalizing, one 
obtains two nearly degenerate neutrino eigenstates with opposite CP parities with a common 



mass of m = ymf + If M is near the GUT scale, m ~ 0.03 eV can be explained quite 
naturally. A bimaximal pattern of neutrino mixings will also be induced, with U given as 
in Eq. (3). 

At this stage, we cannot account for solar neutrino oscillation since Amg = from Eq. 
(6). Once we include the Planck scale corrections from the last term of Eq. (8), one gets a 
mass matrix 



= mi d €2 (9) 




with 5, ej ~ 2 X 10^^ eV. This leads to Am|, ^ 10~^ eV^, which is close to the required value 
for large mixing angle solar neutrino oscillation fit. It is encouraging that all the parameters 
are in the right range to fit solar and atmospheric neutrino oscillation given the assumption 
that the seesaw scale M is near the GUT scale. It should be noted that in this scenario, the 
solar oscillation angle is given by sin^ 26*0 = 1 — O^Am^/ Am^^^Y, which is very close to 1. 
Thus this scenario would prefer the LOW solution to the solar neutrino puzzle, although the 
LMA solution is not entirely excluded. As we shall see in the next section, small breaking 
of Le — Lfj^ — Lr symmetry in the charged lepton sector would make this scenario nicely 
consistent with the LMA solution. 

We now turn to making this type of models more predictive. 



III. RADIATIVELY INDUCED SOLAR NEUTRINO OSCILLATIONS 

In this section, we present a single parameter model for neutrino masses where the 
solar neutrino oscillations are generated by radiative renormalization group effects. We 
will provide a gauge theoretic derivation of the model in section IV. Consider the following 
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neutrino mass matrix defined at a (unification) scale much above the weak scale as in Eq. 
(6). 



Ml= 10 0m. (10) 




Let us assume that the charged lepton mass matrix in the same basis is of the form: 



/ , 

y I . (11) 



\x' 1 

For the most part we will take x' = x so that |x| ~ y^me/m^ and y ~ m^/rrir. We will also 
consider the possibility that x ^ x' (the right handed singlet leptons multiply the matrix 
in Eq. (11) on the right). In the case x = x', note that there is only parameter (m of Eq. 
(10)) in the leptonic sector (both the charged leptons and neutrinos). We will show that 
this model can lead to a realistic description of the neutrino oscillations. 

Below the unification scale, through the renormalization of the effective d = 5 neutrino 
mass operator the form of M° [|T0| will be modified. (Analogous corrections in is 



negligible.) The modified neutrino mass matrix at the weak scale is given by 

- Ml + ^ln{Mu/Mz) {y,Y,^M^^ + M°(r,r/)^) . (12) 

Here c = —3/2 for SM while c = 1 for SUSY, and Yi is the charged lepton Yukawa coupling 
matrix. We have absorbed the flavor-independent renormalization factor into the definition 
of m in Eq. (10). Explicitly, 

2rx 1 + r{x'^ + y'^) 

Mu ^ \ I + r^x"^ + y"^) rx \m, (13) 

rx 



l + r(l + 2x^] 




where 



c 



167r^ 



Y^MMu/Mz) . (14) 



Numerically, r ~ -3.7 x 10"^ in the SM, while it is r ~ 0.022(tan/?/30)2 in SUSY. The 
parameters x ~ 0.017 and y ~ 0.059 are also much smaller than 1, facilitating an expansion 
in small [r, x, y]. Note that in the case of x' ^ x, the numerical value of x can be larger than 
0.017 by a factor even as large as 10. (The e — r mixing will be still small if x ~ 0.17, and 
will not upset the success of approximate Lg — — Lr symmetry.) 
The eigenvalues of to leading order in (r, x, y) are (in units of m) 

2 

mi ^ v^+^(4x + v^) + ^y2 
2 o 

2 

ma ~ -72 - ^(-4x + v^) - ^^2 
2 o 

ms ~ — r^x , (15) 
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where terms of order r^, x^, etc have been dropped. 

The leptoninc mixing matrix U is obtained from U = VeVj , where Vi diagonahzes the 
charged lepton matrix and V^, the neutrino mass matrix: 



U 



i-|(l + v^x) + f| i + i(_l + V2x) 

i + ^ + |(l + v^a:)-^r2 i - ^ + |(1 - v^a:) - ^r^ 
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(16) 



U has approximately the bimaximal form (compare Eq. (16) with Eq. (3)). Let us 
identify mi and m2 to have masses of order 0.05 eV {\/2m ~ O.OSeV) so that the atmospheric 
neutrino oscillation is explained via — oscillations. The solar mass splitting is then 
given by Amg = m\ — m\ ~ 8-\/2rxm^ ~ 0.014rx eV^. In the SM, this is equal to 0.92 x 10~^ 
eV^ if X = x' . On the other hand, if x = lOx', Amg ~ 0.92 x 10"'' eV^, which is in the 
right range for LOW solution of the solar neutrino data. In SUSY, even with x = x', we 
get Amg ~ 4.7 x 10~^ eV^ for tan/? = 30, which is the right range for LMA solution. In 
the case of SUSY with tan/? = 30, if we choose x = 0.2, so that t/gs ~ 0.15, which is near 
its experimental upper limit, we have Amg ~ 5.5 x 10~^ eV^ and the solar mixing angle 
is given by sin^26'0 ~ 0.88. These values are nicely consistent with the large mixing angle 
MSW solution. 

Neutrinoless double beta decay occurs in this model with an amplitude proportional to 
(see Eq. (16)) 



effective ^ 




7712 — \/2xmi ^ 0.05 X eV. (17) 



For X = x' ^ 0.017, this is ~ 10 ^ eV, while for x = lOx', this is 0.01 eV. 



IV. A GAUGE MODEL 

The model for neutrino masses described in the previous section can be derived by 
extending the standard model to include extra gauge singlets 5*0 = (S'l, 5*2) which transform 



as a doublet under an 5*3 permutation group [|Tl|]. The left handed lepton doublets L = 
(L2, L3) of the second and third generation are assumed to transform also as doublets of S3. 
The right handed lepton fields as well as all other standard model fields are assumed to be 
singlets under S^. To obtain the desired charged lepton mass matrix we include two more 
pairs of standard model singlets, Su,d each transforming as a doublet under the 5*3 group. 
We also assume a Z2 symmetry under which 5*^, r^? are odd and all other fields are even. 
The most general Lf. — L^ — Lj. and 5*3 invariant coupling of leptons is given by: 

£1 = -^Le4>L(j)So + -^L(f)SdTR + L(l)Su^i-R + h.c. (18) 

Here M is a fundamental scale close to the GUT scale. The singlets 5*0, Sd all acquire 
VEVs of order M. The 5*3 invariant potential for these fields will admit the following 
structure of VEVs: {Sq) = Vq (1, 1), (Su) = (1, 0) and {Sd) = va (0, 1) [0. Such a VEV 
structure will lead to the charged lepton and neutrino mass matrices given in Eq. (lO)-(ll) 
except for the fact that x = x' = 0. To induce x,x' 7^ 0, we need to introduce singlet Higgs 
fields that break the L^. ~ — symmetry. This can also be done in simple ways. 
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V. A TWO PARAMETER MODEL FROM PERMUTATION SYMMETRY 



In this section, we will present a two parameter model that can also emerge from the 
model of the previous section if we allow for terms that break — — symmetry. 
Consider the following Lagrangian: 

Cl = ■^Le4>L(j)So+ ^{L^,(f)L^(f) + Lr(f)Lr(f)) + L(j)SdTR + L(j)SufJ^R + h.c. (19) 

In this case, one obtains the Majorana neutrino mass matrix to be 

/O 1 1\ 

Ml = \l a m . (20) 
\l a) 

The charged lepton mass matrix is diagonal here. This mass matrix can be diagonalized to 
give the neutrino mixing matrix U, exactly as in Eq. (3). The solar neutrino oscillation is 
governed by AuIq ~ 2m? a. This model predicts t/es = as can be seen from Eq. (3). Thus 
any evidence for nonzero f/es would rule out this model. 



VI. MIRROR NEUTRINOS AND AN EXPLANATION OF LSND EXPERIMENT 



So far we have ignored the LSND results |13] in our discussion of the positive neutrino 



oscillation signals. In order to accommodate the LSND results, one will have to include an 
extra neutrino species that is a standard model singlet (a sterile neutrino). In the presence 
of a sterile neutrino, there are two ways to understand the observations: one is the so- 



called 2+2 scheme 1113 and the other is the 3+1 scheme |T3f. The 2+2 scheme has the 



v^^r neutrinos with mass around an eV and z/e,s with mass near 10~^ eV, with the latter 
explaining the solar neutrino data, the former explaining the atmospheric neutrino data and 
the gap between the two pairs explaining the LSND results. Recent SNO data disfavors the 
original version of the 2+2 model where all the missing solar z/g's are converted via a small 
angle MSW mechanism only to the sterile neutrinos. Mixed 2+2 scenarios where there is 
substantial mixing between the two sectors separated by a gap seem to be consistent with 
all data |TB|. 



In the 3+1 picture, on the other hand, it is assumed that the three active neutrinos are 
bunched together at a small mass value (say around 6 x 10~^ eV or so), with the sterile 
neutrino at a mass near an eV. The atmospheric and solar neutrino data is explained by 
the oscillations among active neutrinos whereas the LSND data is explained by indirect 
oscillations involving the sterile neutrino [1^. The advantage of this picture in view of the 



SNO data is that the flat energy spectrum is explained by postulating bimaximal neutrino 
mixing pattern among the active neutrinos or the LOW solution, as discussed in previous 
sections of this paper. The question now is whether one can extend the model we analyzed 
to include one or more light sterile neutrinos and maintain at least some predictive power 
that we had in the three neutrino case. 

One way to obtain light sterile neutrinos is to invoke a mirror sector for the particles and 
forces, as an exact replica of the standard model sector. One may then identify the neutrinos 
of the mirror sector as sterile neutrinos. The ultralightness of the sterile neutrinos is then 
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explained in a simple way by a mirror version of the seesaw mechanism. In this picture, the 
two sectors are connected only by gravity or possibly by very heavy particles that decouple 
from low energy physics. The mixing between neutrinos then owes its origin either to gravity 
or to operators generated by the exchange of the heavy particles that connect the two sectors. 

In such theories there are two possibilities for the weak scale for the mirror sector: (i) 
it is the same as in the familiar sector or (ii) it is different. In order to explain the LSND 
results in the 3+1 scenario, we need to work with the second scenario, where mirror weak 
scale v' 3> V. We will consider a mirror embedding of the model discussed in section 2 
which was based only on Lf, — L^ — symmetry. We will assume that there is an identical 
symmetry i.e., L'^ — L'^ — L'^ in the mirror sector and as in sec. II we will further assume 
that only Planck scale effects break these two symmetries. They are also responsible for 
mixing between the two sectors. The Lagrangian including the Planck scale breaking terms 
can then be written as: 



+ 
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J a 



al3 



PI 



{L'JL' 



+ bLe(j)Lr(j)) + 



M, 



pe 



IVI " - 



Pi 



+ h.c. 



(21) 



As mentioned before, we work in a picture where the gauge symmetry breaking is not mirror 
symmetric. This then allows for the possibility that both the weak scale and the seesaw scale 
in the mirror sector are very different from those in the familiar sector. Suppose we assume 
that {(/)') = i;' ^ 5 X 10^ GeV and the mirror seesaw scale M' ^ 10^ GeV or so. From 
Eq. (21) one sees that the two heavy neutrino states z/( 2 can have masses easily of the 
order of 1-10 GeV and the lightest mirror neutrino has a mass of ~ eV and can therefore 
be used to understand the LSND results. The off diagonal v — v' mixing terms are of order 
vv' jMpi ~ 0.01 eV, which is also of the right order for the purpose. 

Upon integrating out the heavy neutrinos, one gets a 4x4 matrix with the following 
generic structure (neglecting terms smaller than 10~^ eV): 
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mi 


m2 




mi 








£2 


m2 








£3 


V ei 


£2 


£3 


5/ 



(22) 



For b ^ €i,mi, this matrix has one large eigenvalue ~ 6, which mixes with u^^^ with a mixing 
angle of order ei/6; for ~ 0.02 eV and S ~ 0.2 eV, the LSND results can be explained 



using indirect oscillation [|T8[, with the probability given by P(z/e 



1^1, 



4|eie2/5 



212 



VII. COSMOLOGICAL IMPLICATIONS 

The mirror sector scenario has several interesting cosmological implication. First of all, 
the two heavy GeV) mirror neutrinos annihilate into the lightest mirror neutrino early 
on and disappear. As a result they do not contribute to the expansion of the universe 
during big bang nucleosynthesis (BBN). The lightest mirror neutrino will however be in 
equilibrium during BBN. Thus, this model would predict one extra neutrino species at the 
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BBN. In principle, this could be tested once the He^ and abundances are known more 
precisely. 

A second interesting implication is that due to the heavy mirror scale, the mirror QCD 
scale is A' !^ 7 GeV (for the nonsupersymmetric case). This is calculated as follows. Mirror 
symmetry guarantees that the two QCD couplings are the same at very high energies. As 
they evolve to lower scales, due to the heavier masses of the mirror quarks a'^ will develop 
a higher slope around 10^ GeV corresponding to t' decoupling and become larger than a^; 
this effect gets further reinforced at the (6', c') masses. A simple calculation using as{Mz) — 
0.118 shows that a'^ ~ 1 around 7 GeV. We define this as the mirror QCD scale. Thus from 
the neutrino oscillations we are able to predict the mirror QCD scale. Using the analogy 
with familiar QCD then one can predict that mirror proton mass is around 30-40 GeV. 

The mirror electron mass is mg' ~ — 500 MeV. This would make the mirror 
hydrogen Bohr radius 1000 times smaller than the familiar hydrogen atom, leading to a 
completely new profile for the mirror sector of the universe. 

The mirror baryons, being stable, can now become the dark matter of the universe. This 



picture is very similar to the one advocated in Ref. |]T9|, with an important difference that 



here the mirror electron is about 10-30 times heavier. To see how this fits in with the dark 



matter picture pO[, we note that in this scenario, to reconcile the presence of the mirror 



photon with the success of BBN, one needs to have the temperature of the mirror universe 



somewhat (say a factor of 2) lower than the temperature of the familiar sector [^. One 
can then write 



(23) 



For nipt /nip ^ 40 and T' /T ~ 0.5, we get ^Ib'/^b — 5. If we take Hb — 0.05, this leads to 
floM = flB' ^ 25%, consistent with the current thinking in the field of dark matter physics. 



VIII. CONCLUSIONS 

We have presented several simple and predictive three neutrino oscillation scenarios that 
can explain the solar and the atmospheric neutrino data with a bimaximal mixing pattern 
that seems to be favored by current data, especially from the fiat energy distribution for 
solar neutrinos. The measurement of the mixing parameter Ues and evidence for inverted 
pattern of masses would constitute tests of these models. We have also presented a mirror 
extension of the schemes that can accommodate the LSND results as well. We find that 
in the simplest mirror extension, the mirror QCD scale is predicted to be around 5-7 GeV. 
This leads to a mirror baryon mass mass around 40 GeV, which can then become a viable 
candidate for dark matter of the universe. 
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